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A BCS Gap on the Lattice 

David N. Walters^ and Simon Hands^* 

'^Department of Physics, University of Wales Swansea, 
Singleton Park, Swansea SA2 8PP, United Kingdom. 

Monte Carlo simulations of the 3+1 dimensional NJL model are performed with baryon chemical potential 
jj, > 0. For /i ^ Eo, the constituent quark mass in vacuum, chiral symmetry is restored and a diquark condensate 
{qq+) forms. We analyse the fermion propagator and find evidence for particle-hole mixing in the vicinity of 
the Fermi surface and an energy gap A > 0, both of which provide evidence for superfluidity at high baryon 
density induced by a BCS mechanism. At fia — 0.8 the ratio between the BCS gap and the vacuum quark mass 
is A/Eo = 0.15(2). 



The Nambu - Jona-Lasinio (NJL) model has 
served for many years as a model of strong in- 
teractions PI; in particular with the addition of 
a non-zero baryon chemical potential /i it can be 
used to model the relativistic quark matter an- 
ticipated in dense systems such as neutron star 
cores. Phenomena such as the BCS-condensation 
of quark pairs at the Fermi surface, suggested 
as a mechanism for color superconductivity and 
baryon number superfluidity, can be modelled in 
this way 0. From a lattice theorist's perspec- 
tive, the attractive feature of the NJL model is 
that its action in Euclidean metric remains real 
even once ^ ^ 0, so that standard Monte Carlo 
simulation techniques can be employed. More- 
over, such simulations expose a chiral symmetry 
restoration transition at a critical /^c ~ So, where 
Eo is the constituent quark scale, in qualitative 
agreement with analytic approaches such as the 
Hartree approximation 

In previous work [3] we have simulated the NJL 
model in an attempt to identify diquark conden- 
sation for jj, > ^c- In continuum notation the 
Lagrangian is 



C = + m + Ai7o)V' 



(1) 



,tr\ I J 



J J V V' 



When formulated using staggered fermions, the 
resulting lattice model has Nf = 2 isospin degrees 
of freedom, and Nc = 8 global 'colors' (there are 
no gluons). The third term in is a U(1)b- 
violating source term which enables the identifi- 
cation of a superfluid diquark condensate {qq+) 
for source strengths j — J = j* 7^ on a finite 
volume. 

At /i = the model exhibits breaking of 
SU(2)l«)SU(2)/j into SU(2)/ by the spontaneous 
generation of a chiral condensate and a con- 
stituent quark mass Sq ^ m. Since in 3+1 
dimensions, the NJL model is only an effective 
theory, the lattice parameters must be chosen 
to match low energy phenomenology. Using the 
large- iVc (Hartree) approximation we determine a 
suitable set, leading to physically reasonable re- 
sults, to be 
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ma = 0.006 
2 = 0.495 
1 = 720MeV 



Eo = 400MeV 
U = 93MeV (2) 
= 138MeV. 



We use a standard HMC algorithm with a quark 
kinetic matrix M which is a functional of auxil- 
iary boson fields a and 7?. The source strength j is 
set to zero in the HMC update, but is allowed to 
vary over a range of values in partially quenched 
measurements on valence quarks. If we define an 
antisymmetric matrix 

1 / 7T9 M 
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(3) 
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then e.g. the diquark condensate is given by 

, . 1 9 In Z 1 , , , , , ^ , 

(qq,) = — = trra^"^). (4) 

The results show that at ficO- — 0.6 there is 
a transition at which the chiral condensate de- 
creases sharply until by = 0.8 it has just 5% 
of its vacuum value, and at the same point the 
baryon density ub = {2V)^^d\n Z /dfi begins to 
rise from zero. Results for {qq+), first extrap- 
olated linearly in Lj"^ to T = from 12'', 16"* 
and 20^ lattices, and then linearly in j — > using 
data with ja £ (0.3,1.0), suggest that a super- 
fluid condensate forms for fia ^ 0.6 and then in- 
creases monotonically until /xa ~ 1.0, where sat- 
uration artifacts become apparent 0. 

In a system with a Fermi surface, low-energy 
excitations have k = \k\ close to the Fermi mo- 
mentum hp. If a BCS condensate forms, we ex- 
pect an energy gap 2A to open up between the 
highest occupied state in the Fermi sea and the 
lowest excited state. These excitations are probed 
by standard means via temporal decay of the Eu- 
clidean propagator Q = A'^. We write 

G(x v)- f ^^""'^^ ^i^^y)\ (.) 
^^""^y) ' [N{x,y) A{x,y))^ 

where we identify both "normal" N{x,y) ^ 
{q{x)q{y)) and "anomalous" A{x,y) ^ {q{x)q{y)) 
propagators. In an isospin-symmetric ground 
state the only independent components to sur- 
vive the quantum average are N = ReiVn and 
A = ImAi2, where subscripts label isospin com- 
ponents We proceed by studying the time- 
slice propagator g{k,t) = ^^(0, 0; f , i)e-'^ =^ 
on a 96 X 12^ x Lt lattice {Lt = 12, 16, 20 and 
24) with k = (7rn/48, 0, 0), n = 0, 1, . . . , 24. The 
resulting propagators are fitted with 

N{k,t) = Ae-^* -f- Be--^(-^'-*) 

A{k,t) = C(e-^* -e-^(^'-*)) (6) 

and amplitudes A, _B, C, and excitation energy 
E{k,j) extracted. 

Fig. n shows the pole-fit amplitudes as func- 
tions of k. For small k the normal propagator 
is predominantly forward- moving, since A B. 
For k ^ fj,, however, the signal becomes predom- 
inantly backward-moving. This is a sign that the 
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Figure 1. Amplitudes at fia = 0.8 extrapolated 
to T followed by j ^ 0. 



excitations for small k arc hole-like, and for large 
k particle- like 4 . Of particular interest, however, 
is that in this case even in the limit j there 
is a region of fc « /i where the anomalous ampli- 
tude C ^ 0. This shows that there are propagat- 
ing states which are particle-hole superpositions 
with indefinite baryon number, an indirect signal 
for superfluidity via a BCS mechanism. 

In Fig. [5] we show the dispersion relation E{k) 
extracted from the anomalous propagator at fxa = 
0.8. It is important to note the order of extrap- 
olation: first r ^ 0, then j 0. For clarity 
we have plotted energies for fc < kp as negative, 
where the Fermi momentum fci? is chosen to be 
the point where ^ = i? in Fig. ^ We have also 
plotted the free lattice fermion dispersion rela- 
tion, which has i?(sin~^(sinh(/ia)) = 0. The in- 
teracting data, however, show no sign of a zero 
energy state; instead there is a a discontinuity at 
k K, kp signalling an energy gap 2 A. This is the 
first direct observation of a BCS gap in a lattice 
simulation. 

Finally in Fig.|31we plot both order parameter 
{qq+) and gap A as functions of /i. Note that a re- 
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Figure 2. Dispersion relation E{k) at /la = 0.8 



Figure 3. {qq+) and A as functions of /i. 



liable extrapolation to T — > requiring data from 
Lt = 24 lattices is only available for jia = 0.8; 
simple linear extrapolation through Lt = 16 and 
20 is completely consistent with this, however, 
and hence is used for the other ^ values. One 
physical interpretation of the order parameter is 
that it counts the number of qq pairs participat- 
ing in the condensate. On the simple-minded as- 
sumption that these pairs occupy a shell of thick- 
ness ©(A) around the Fermi surface, we surmise 
the relation 



{qq+) oc A^^ 



(7) 



Fig. 121 indeed suggests that A is roughly constant 
for fi > fic, while {qq+) rises markedly. The ratio 
A/Eq ~ 0.15, which using Q translates into a 
physical prediction A « 60MeV, consistent with 
the model predictions of j^]. 

In summary, by presenting numerical estimates 
for both {qq+) and A we have amassed a reason- 
able body of evidence for superfluidity via a BCS 
mechanism in a relativistic quantum field theory, 
for the first time using a systematic calculational 
technique. In future work we hope to address the 
issue of finite volume effects, unusually large in 
this system which we believe is due to the 
difficulty of representing a curved shell of states 



around the Fermi surface on a discrete momen- 
tum lattice. It will also be interesting to study 
the stabihty of the superfiuid phase for T > 0, 
and for /x„ ^ /i^ , which has the effect of separat- 
ing the Fermi surfaces of each flavor and hence 
suppressing isosinglet pairing. 
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